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Abstract. To obtain an exact solution of a four-body system containing two 
quarks and two antiquarks interacting tiirougli two-body terms is a cumber- 
some task that has been tackled with more or less success during the last 
decades. We present an exact method for the study of four-quark systems 
based on the hyper spherical harmonics formalism that allows us to solve it 
without resorting to further approximations, like for instance the existence of 
diquark components. We apply it to systems containing two heavy and two 
light quarks using different quark-quark potentials. While QQnn states may 
be stable in nature, the stability of QQnn states would imply the existence of 
quark correlations not taken into account by simple quark dynamical models. 

The discoveries on several fronts [Ij, of unusual charmonium states like 
X(3872) and y(4260) and open-charm mesons with unexpected masses like 
L'*j(2317) and Dq{2?>QS), have re-invigorated the study of hadronic resonances. 
Any debate on the possible multiquark structure of meson resonances should be 
based on our capability to find an exact solution of the four-body problem [2]. 
Theoretical predictions often differ because of the approximation method used. 
A powerful tool to solve a few-particle system is an expansion of the trial wave 
function in terms of hyperspherical harmonics (HH) basis functions. In Ref. [3] 
a generalization of the HH formalism to study four-quark systems in an exact 
way was presented. Due to their actual interest and having in mind that sys- 
tems with unequal masses are more promising [2], we will center our attention 
on QQnn and QQnn states (n stands for a light quark and Q for a heavy one). 
We will analyze the possible existence of compact four-quark bound states using 
two standard quark-quark interactions, a Bhaduri-like potential (BCN) [3j and 
a constituent quark model considering boson exchanges (CQC) [5j. Both inter- 
actions fulfill the requirement of giving a reasonable description of meson and 
baryon spectroscopy. Assuming non-relativistic quantum mechanics we solve the 
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four-body Schrodinger equation. The grand angular momentum K is the main 
quantum number in our expansion and the calculation is truncated at some K 
value. Further details of the numerical method can be found in Ref. [3]. 

In spite of the shortcomings of the methods used to study four-quark systems, 
in the past, many four-quark bound states have been suggested. To analyze their 
stability against dissociation, parity and total angular momentum must be pre- 
served. Additionally, C— parity is a good quantum number for ccnn and the Pauli 
principle must be fulfilled in the ccnn case. The corresponding thresholds can be 
computed by adding the meson masses of the dissociation channel. Four-quark 
states will be stable under strong interaction, and therefore very narrow, if their 
total energy lies below all allowed two-meson thresholds. Sometimes, results of 
four-quark calculations have been directly compared to experimental thresholds. 
In this case one could misidentify scattering wave functions as bound states. 
When they are referred to the thresholds within the same model, theoretical 
predictions do not imply an abundance of multiquark states in the data. 

Table 1. Energy (MeV) and probability of the different color components 
for the ccnn J^'^ = 1++ both for CQC and BCN models. The last rows 
indicate the lowest theoretical two-meson thresholds. Pn (Pss) stands for 
the probability of singlet-singlet (octet-octet) color components. 
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Once the method has been established, we concentrate on the ccnn systems 
as a potential structure for the X(3872). To make the physics clear we com- 
pare with the ccnn system. In particular, we focus on the J^'^ = l"*"*" ccnn 
and = 1+ ccnn quantum numbers to illustrate their similitudes and dif- 
ferences. A complete study of all the quantum numbers have been reported in 
Ref. [3J . For the ccnn system, independently of the quark-quark interaction, the 
system evolves to a well separated two-meson state, see Table . This is clearly 
seen in the energy, approaching the corresponding two free-meson threshold, 
but also in the probabilities of the different color components of the wave func- 
tion. Comparing the theoretical predictions with the experimental threshold, 
^J/tpLu\s ~ 3879.57 ±0.13 MeV, one could be tempted to claim for the existence 
of a bound state. However, the experimental threshold is not reproduced by the 
effective Hamiltonians. Thus, in any manner one can claim for the existence of 
a bound state. Similar conclusions are drawn for all quantum numbers of this 
system. A completely different behavior is observed in the case of = 1"*" ccnn. 
The energy quickly stabilizes below the lowest theoretical thresholds (3937 MeV 
for CQC and 3906 for BCN), being the results obtained for Kmax = 24 com- 
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pletely converged, 3861 MeV for CQC and 3900 for BCN. Besides, the radius is 
also stable and is smaller than the sum of the radius of the two-meson threshold. 
We obtain r^g = 0.37 fm compared to tmi + ^"Ah = 0-44 fm. 

It is thus important to realize that a bound state should be pursued not 
only by looking at the energy, but also with a careful analysis of the radius and 
probabilities. This detailed analysis allows us to distinguish between compact 
states and meson-meson molecules [6] and it does consider the contribution of 
all meson-meson channels to a particular set J^*^ of quantum numbers. Inherent 
to our discussion is a much richer decay spectrum of compact states due to the 
presence of octet-octet color components in their wave function. 

Let us notice that there is an important difference between the two phys- 
ical systems studied. While for the ccnn there are two allowed physical decay 
channels, {cc){nn) and (cn)(cn), for the ccnn only one physical system contains 
the possible final states, (cn)(cn). This has important consequences if both sys- 
tems (two- and four-quark states) are to be described within the same two-body 
Hamiltonian, the ccnn will hardly present bound states, because the system will 
reorder itself to become the lightest two-meson state, either {cc){nn) or {cn){cn). 
In other words, if the attraction is provided by the interaction between parti- 
cles i and j, it does also exist in the asymptotic two meson state reflecting this 
attraction. This may not happen for the ccnn if the interaction between, for ex- 
ample, the two quarks is strongly attractive. In this case there is no asymptotic 
two-meson state with such attraction, and therefore the system will bind. 

Therefore, our conclusions can be made more general. If we have an A^-quark 
system described by two-body interactions in such a way that there exists a 
subset of quarks that cannot make up a physical subsystem, then one may expect 
the existence of A^-quark bound states by means of central two-body potentials. 
If this is not true one will hardly find A^— quark bound states [7j. 
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